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Abstract 

Analytic formulae are given for order a s perturbative QCD correc¬ 
tions to the double differential distribution of the r lepton energy and 
the invariant mass of the r + antineutrino system in semileptonic de¬ 
cays of the bottom quark. The corresponding distribution for B meson 
decays are obtained by combining these results and the nonperturbative 
QCD corrections calculated by Falk et al. The moments of the energy 
distribution of the r lepton are calculated. 
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1 Introduction 


The theoretical description of the semileptonic decays of beautiful hadrons 
based on the Heavy Quark Effective Theory (HQET) |], ||, ff 0. |], § or 
the heavy quark mass expansion ]7|, ||, |Dj| and on perturbative QCD seems 
to be quite satisfactory. The recent ALEPH measurement of the semilep¬ 
tonic branching ratio for the baryon Afe|10[ indicates that the problems of the 
present theory to explain the short lifetime of A;, are associated with the sector 
of nonlcptonic decays. In contrast to the nonleptonic widths, the semileptonic 
widths of beautiful hadrons seem to follow the theoretical predictions. This 
conclusion is strengthen by an old observation that, despite the large differ¬ 
ences in the lifetimes, the semileptonic widths of charmed hadrons are fairly 
similar. It is clear that the theory of the hadronic processes involving heavy 
quarks is still to be improved. Another failure is the prediction of a large A& 
polarization in Z° decays. In fact the measurement HTTfl of A^ polarization via 
its semileptonic decay channels |[T2|. [Pij . |Tj|j gives a much smaller value than 
predicted in the heavy quark mass limit ||15||. All this indicates that further 
developments in the theory of nonleptonic processes are expected whereas the 
theory of semileptonic decays is firm and stable. Therefore the latter pro¬ 
cesses can be used in determination of parameters of the standard model like 
the masses mb and m c of the charm and bottom quarks and the strong cou¬ 
pling constant a s at relatively low scales. Recent determination [|l6j of the 
Cabibbo-Kobayashi-Maskawa matrix element V c b is a spectacular example. It 
has been argued [|13|, [L7|, |l|] that moments of the lepton spectra can be used in 
determination of a s , nib and m c . Semileptonic decays of the bottom quark are 
particularly promising in this respect because in addition to the channels with 
practically massless charged leptons (electron and muon) in the final states 
the channel is kinematically allowed with the r lepton, whose mass is appre¬ 
ciable in comparison to mb and m c . This process can be also a window on 
a new physics Jl^, 


21], The theoretical description of inclusive semilep¬ 


tonic decays of the b quark into the r lepton in the framework of the standard 
model is therefore necessary. The nonperturbative QCD corrections to the 
total rate and the r lepton energy spectrum have been calculated by a number 
of groups[^2], ^3], ^4|. Perturbative QCD corrections to the total semileptonic 
decay rate are known. In |22| a numerical approach was used based on an 
important work by Hokim and Phamp5|; see also Bagan et al 126[. In addi¬ 


tion in|17[ an analytical formula has been obtained for the distribution of the 


1 








































invariant mass of the virtual W (i.e. the tD t system) which after performing 
one-dimensional numerical integration gives the correction to the semilcptonic 
rate in perfect agreement with the results of [f£^]. The case of perturbative 
QCD corrections to the r lepton energy spectrum is less satisfactory. In ref. |27j 


formulae are given for these corrections as a one-dimensional integral. In the 
present article we give another representation of these corrections. We have 
calculated analytically the double differential distribution of the r energy and 
the mass of the tD t system. In the limit m T —> 0 our formulae reduce to the 
well-established results for the massless case p8|, |29| . After numerical integra¬ 
tion over the energy of r they are in perfect agreement with the analytical 
formula given in[[L7]]. We were unable, however, to find agreement between our 
results and those of ref. [[jF/J for the corrections to the r energy distribution. 
Our numerical results for the r energy distribution have been given in |3(J. In 
the present article we calculate the moments of this distribution. 

The article is organized as follows: in Sec.2 kinematical variables are de¬ 
fined which are used throughout. In Sec.3 the calculation of QCD corrections 
is briefly described. In Sec.4 we give the analytic result for the double differen¬ 
tial distribution of the r energy and the invariant mass of the tv t system. This 
is the main result of the paper. In Sec.5 the moments of the r lepton energy 
distribution are calculated numerically and tabulated. In Sec.6 the results of 
this article are summarized. 


2 Kinematics 

2.1 Kinematical variables 

The purpose of this section is to define the kinematical variables which are used 
in this paper. We describe also the constraints imposed on these variables for 
three and four-body decays of the heavy quark. 

The calculation is performed in the rest frame of the decaying b quark. 
Since the first order perturbative QCD corrections to the inclusive process are 
taken into account, the final state can consist either of a produced quark c, 
a lepton r and an anti-neutrino v T or of the three particles and a real gluon. 
The four-momenta of the particles are denoted in the following way: Q for the 
b quark, q for the c quark, r for the charged lepton, v for the corresponding 
anti-neutrino and G for the real gluon. By the assumption that all the particles 
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are on-shell, the squares of their four-momenta are equal to the squares of 
masses: 


Q 2 = ml 


Q 2 = ml 


2 2 
r = m T , 


V 2 = G 2 = 0. 


( 1 ) 


The four-vectors P = q + G and W = r + v characterize the quark-gluon 
system and a virtual intermediating W boson respectively. We define a set of 
variables scaled in the units of mass of the heavy quark mb'- 


P = 


mi 


m 


2 ’ 


V = 


mz 


2 E t 


m: 


2 ’ 


X = 


t = 


W 2 


P 2 


mi 


m 


2 ’ 


2 = 


mi 


( 2 ) 


We choose as the unit of mass (i.e. Q 2 
variables describing the charged lepton: 



1) and introduce light-cone 


( 3 ) 


The system of the c quark and real gluon is characterized by the following 
quantities: 


PoM 

p 3 M 
p±W 


-(1 — t + z), 

+i 2 + 2 2 


P 0 (*)±^(*), 

i ln ^44 =ln ^M, 

2 P_(z) ^ ’ 


2(f + ^ + f^)] 1 / 2 , 


( 4 ) 


where Po(z) and P 3 {z) are the energy and the length of the momentum vector 
of the system in the b quark rest frame, y p (z) is the corresponding rapidity. 
Similarly for the virtual boson W: 


W 0 (z) 

W 3 (z) 

W±(z) 

y w (z ) 


2 (1 + t ~ z )> 

\JWq — t = - [1 + t 2 + z 2 


W 0 (z)±W 3 (z), 

(z) 


iln^- 

2 W-(z) 


= In 


W+(z) 

Vi 


2 (t + z + tz)} 1/2 , 


( 5 ) 
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^Frorn kinematical point of view the three body decay is a special case of the 
four body one with the vanishing gluon four-momentum, what is equivalent to 
a substitution z = p in the previous formulae. It is convenient to use in this 
case the following variables: 

Po = Po(p) = ^(1 ~t + p), Pi = P?Ap) = yj'pi - p, 

P± = P±{P) =P0± P3, W± = W±(p) = 1 - p -. F , 

U = %(p) = bn—, U = y»(p) = bii—■ (6) 

2 p_ 2 

We express also the scalar products which appear in the calculation by the 
variables x, t and z: 

QP = |(1 + z-t), t-v 

Q-v — |(1 - z - X + t), T-P — \{x - t - Tj), (7) 

Q-r — \x, v-t — |(1 — x — z + rj). 

All of the written above products are scaled in the units of the mass of the 
b quark. 

2.2 Kinematical boundaries 

The allowed ranges of x and t for the three-body decay are given by inequali¬ 
ties: 


2 y/rj <x<l + rj-p = x max , (8) 

il = T _( 1 - r ^-)<i< T + ( 1 - r iL_)= i2 (9) 

(a region A). In the case of the four-body process the available region of the 
phase space is larger than the region A. The additional, specific for the four 
body decay area of the phase space is denoted as a region B. Its boundaries 
are the following: 

2y/p <x< x max , rj <t < ti. (10) 
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We remark, that if the charged lepton mass tends to zero than the region B 
vanishes. 

One can also parameterize the kinematical boundaries of x as functions 
of t. In this case we obtain for the region A: 

V < t < (1 - v/p) 2 , w- + — < x < w + + —, (11) 

w _ w + 

and for the region B: 

V < t < Vv ^ x - ’ 2 ^ Tl ~ x ~ w ~ + ( 12 ) 

The upper limit of the mass squared of the c-quark — gluon system is in 
the both regions given by 

Anax = (1 — t+)(! ~t/r + ), (13) 


whereas the lower limit depends on a region: 


p in the region A, 

(1 — r_)( 1 — t/rJ) in the region B. 


(14) 


3 Calculation of QCD corrections 

The QCD corrected differential rate for b —> c + r~ + V reads: 

dT = dT 0 T dT 13 T dTi 4 , (15) 

where 

dr 0 = G 2 F ml\V C KM\ 2 Mo t 3 dn 3 (Q-, q, r, z/)/r 5 (16) 

in Born approximation, 

dTi )3 = ^a s Gprnl\V C KM\ 2 Mi j3 dn 3 (Q-,q,r,u)/-K 6 (17) 

comes from the interference between the virtual gluon contribution and Born 
amplitudes, and 

dr i )4 = ^a s G 2 F ml\VcKM\ 2 Mi j4 d7l4(Q; q, r, zz)/7r 7 (18) 
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describes a real gluon emission. Vckm is the Cabibbo-Kobayashi-Maskawa 
matrix element associated the b to c or u quark weak transition. Lorentz 
invariant n-body phase space is defined as 

dTZ n (P;Pi, ■ ■ ■ ,Pn) = ( 19 ) 

l L 

In Born approximation the rate for the decay into three body final state is 
proportional to the expression 

= F 0 {x,t) = 4q-r Q-u — (1 - p - x + t)(x — t — 77 ), (20) 

where the quantities describing the W boson propagator are neglected. The 
three-body phase space is parameterized by the use of Dalitz variables: 

7r 2 

dH 3 (Q', q, t, v) = —dxdt. (21) 

The matrix element M .13 was evaluated in the fifties |3l[ and has the following 
form: 


M l3 = —[H 0 q-T Q-v + H + p Q-v Q-r + H_ q-v q-r 

+\p( H + + H -) v-t + ^pp(H + - H_ + H l ) Q-v 

-\ H Lpq-v}i 


( 22 ) 


where 


H 0 = 4(1 - Yppo/po ) In A G + (2p 0 /p 3 ) 


H± 

H l 


Lio I 1 — 


\ 

p+w+J 


Lio 


W- 


1 - — I - Y P (Y P + 1) + 2(In yfp + Y p ) (Y w + Y p 


+ [2 p 3 Y p +(l-p-2t) In y/p\/t + 4, 

\\ l ± (1 - P)/t]Y p /p 3 ± j In y/p, 

1(1 - In y/p) + Yl In Vp + ly > 3 + jL. 


(23) 


The virtual correction in the interference term is renormalized and hence Adi )3 
is ultraviolet convergent. However, the infrared divergences are left. They are 
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regularized by a small mass of gluon denoted by A a- According to Kinoshita- 
Lee-Naunberg theorem, the infrared divergent part should cancel with the 
infrared contribution of the four-body decay amplitude integrated over suitable 
part of the phase space. 

The rate from real gluon emission is proportional to 


KA — 1 _=_|_ £ 

1)4 (Q-G ) 2 Q-GP-G (P-G) 2 ’ 

where 


(24) 


B, = q -T [Q-v (Q-G - 1) + G-v - Q-v Q-G + G-v Q-G], 

— q-T [G-v Q-q — q-v Q-G + Q-v (q-G — Q-G — 2q-Q)] 

+ Q-v (Q-t q-G - G-t q-Q), 

— Q-v (G-t q-G — p t-P). 

(25) 


The four-body phase space is decomposed as follows: 

dTZ 4 (Q] q, G, r, v) = dz dTZ 3 (Q] P, r, v) d1Z 2 (P] q, G). 


(26) 


The four-momentum of the c quark is substituted by P — G and the integration 
of -Mf ,4 over dlZ 2 (P\ q,G) is performed. Lorentz invariance allows to reduce 
all of appearing integrals to scalar ones: 

I n = f dn 2 (P-q,G)(Q-G) n . (27) 

The formulae for I n have been presented in p9| . 

In the next step Dalitz parametrization (0) of the three-body phase space 
<777-3 (Q ; P, t, v) is employed. The expression obtained after integrations has 
one part 


const F 0 (x,t)I div , 


(28) 


where 


/div-/ 2 (1 t P)I-i/(P'G) + pIo/(P-G) 2 , (29) 

which is infrared divergent when (x, t) belongs to the region A and the rest 
which is infrared finite. Since the method used in this calculation is exactly 
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the same as this used in the previous ones |28], ^J], we do not go into the 
details. The important point is, that the infrared divergent part is regularized 
by a small gluon mass A q and then contains terms proportional to In Xq- The 
infrared divergent terms from three- and four-body contributions to the decay 
rate cancel out and the limit of vanishing gluon mass Xq —> 0 is performed. 
This procedure yields well defined double differential distributions of lepton 
spectra which are described below. 


4 Analytic results 

The double differential unpolarized distribution of the r energy and tv invari¬ 
ant mass squared from the b quark decay with first order perturbative QCD 
corrections reads [] 


dT 

dx dt 


12T n 


F 0 (x, t ) - t)J for (x, t ) in A, 

for (x, t ) in B, 


12T 0 ^Ff(x,t) 


where 


T n = 


Gpinf 


\V C 


ICM 


and 


192vr 3 

F 0 (x,t ) = (1 - p - x + t)(x - t — rf) 

F^x, t) = F 0 % + D*$ n + Dt 

n=1 


(30) 

(31) 

(32) 

(33) 


Ff(x, t) = f „* 0 + y r>y n + D 


n= 1 


(34) 


The factor of 12 in the formula (]3D|) is introduced to meet widely used |T7j. [22|, 
18| convention for F 0 (x) and T 0 . The symbols present in fl33U34[) are defined 


1 A Fortran version of the formulae given in this section is available upon request from 
jezabek@hpjmiady.ifj.edu.pl or leszekm@thpl.if.uj.edu.pl. 
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as follows 


$0 


$1 

$2 

$3 

$4 

$5 


2 p 0 

P3 


Lio 


I. i; 


1 — r H 
P+ 


P- 


+ Li 2 


Lio 


1 — t/r + 
P+ , 

1 zlhl\ 

P- ) 


+Li 2 (w-) - Li 2 (w+) + 4 Y p In yfp\ 

+4 (1 — —Y p ) In ( 

^max p) — 4 In ^ max , 

\ P3 J 

Li 2 («;_) + Li 2 (w+) - Li 2 (r+) - Li 2 (t/r+), 
Pa’ 

2 ln a/p> 

i ln(l — r + ), 


In (1 -t/r+), 


(35) 




4/ 2 

^3 


Po 

\P3 

2po 

P3 


Y p - 1 ln 


P 


P, 


Li 2 ^1 - 
—Li 2 ^1 - 
+Li 2 (1 — 


1 ~ t+ 

P- 

1 ~T+ 

P+ . 

1 — r_' 


+ 41n(^) 

\ Z m [ n / 

1 -t/T+' 


+ Li 2 1 - 


-Lio 1 


P+ y 
1 — 

P- , 


Li 2 1 


+ Li 2 1 - 


Lio 1 


P- 

1 - */t+ N 

P+ y 
1 — t/r_ N 
P+ y 
1 — t/r_ 
P- 


Li 2 (r + ) + Li 2 (t/r+) - Li 2 (r_) - Li 2 (t/r_), 
Ln(l-r_), 

Ln(l -t/rJ), 


5 
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(36) 


^4 = ^ ln(l — t + ), 
^5 = ^ ln(l — t/r + ). 


) B : 

n 


We introduce C\ ... C 5 to simplify the formulae for D// and D* 

—2 xt + x + t + t 2 — px + pt + 5p + p(x — t — p), 
\Jx 2 — 4p (6 — 2 x + 2 t — t 2 /p + p — 2 p), 

—5 — x + 2a; 2 — 3 t + 5 1 2 — 6 xt + xt 2 /p — 2 px + pt 
+llp — p 2 + p( 4 + 7x — 11 t — 7 p) + p 2 , 

—4 xt — 3x + 2a; 2 + lit + 2f 2 — xt 2 Ip + 2 xp/t — 2 px 
—10 p/t + 3 pt + 3 p + 5 p 2 /t 2 — Qp 2 /t + p(5x — 71 
+ 2 xp/t + 8 p/t — 11 p — Ap 2 /t 2 ) + p 2 (2 — p/t)p/t , 


<?! = 
C*2 = 

^3 = 

c 4 = 


c, = 


■j- 

-\/x 2 — 4 p 5 — 2t 2 /p — 5p/t — p + ‘M — 3p(l — p/t) 
, pp(l + p)(l- pA) p(t- p) 


t x — t — p/t 1 — x + p 


Dt = C u 


D? = 


(—5 + 2x — 3t + p)pl — 2 pt + 2p{—4p\/t 2 + [x — 7t 
+t 2 - p(x + t)\pl/t 2 } + p 2 [(3 + t + p)p\/t 2 + 2p/t], 


£>3 = —4p^ — 4 — 2xt + 2x + 4f 2 + 2p(2 + 3a; — 7t) 


Dt = 

D/ = 


+2//[4(a; — 1 + p)p 2 z /t 2 + 6 + x/t — x — 6 jt 

+p(—6 + x/t + 6 /t)\ + 2// 2 [—2p 2 + 2 — 2 1 — p(2 + t)]/t 2 

-C 2 - C 3 , 

C 2 - c 4 , 


(37) 


Dq = -(5a; + 9xt + 4t — 6f 2 — 2xt 2 /p — pxjt + px — Ap/t — 2pt 


-22p + 6 p 2 /t) + -p(—3a; + 5f — xp/t + 9p/f + 5p 

+v/t) - Ar' )( r ,,) - t ! a< 3 +•»/*) 

4 1 — a; + p 4 t 

1 pp (1 + p)(t - p/t)(l - p/t) _ 1 
4 t x — t — p/t 2 5 ’ 


(38) 
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Df = C u 
D* = C 2 -C 3 , 
Di = -C 2 -C 4 , 
Df = C 2 + C 3 , 
Di = ~C 2 + C' 4 , 
Di = C 5 . 


( 39 ) 


One can perform the limit p —> 0, what corresponds to the decay of the 
bottom quark to an up quark and leptons. The formulae, which are much 
simpler in this case, are presented in the same manner as previously the full 
results: 


dT 

dx dt 


where 

Fo(x,t) = 

and 


12 T 
12 T 


o F a (x,t) - 


2 a s 

3tv 




(1 — x + t)(x — t — rf) 


for (x, t ) in A, 
for (x, t ) in B, 


Fi(x, t ) = F 0 $0 + Ztf *1 + v£ 3 $ 2<>3 + rtf $4 + $5 + Ztf, 


(40) 


(41) 

(42) 


F 1 B (x,t) = FoT 0 + ^P,fT n + D ( 


6 ) 


n=l 


where 


$o 


4>i 


$203 

$4 

$5 


Lb 


T + ~t 


+ Lh 


A/ r +~ 

V !A -1, 


+ Li 2 (i) 


+ 2 ^ 


+ ln 2 (l - r+) + 2 ln 2 (l - f) + ln 2 (l - t/r+) 


—2 ln(l — t) lnz max , 

7T 2 

— + Li 2 (t) - Li 2 (r+) - Li 2 (t/r+), 
2 ln(l — t) 


$4, 

$5 


(43) 


( 44 ) 


11 



and 


'Ln = 2 


Li 9 


r_ — t 
1 -t 

-Li 2 ^ 1/r+_ 


“1“ Li 


l/r_ — l' 


— Lio 


In 

In 


V l/t-l 

'&) 

1 - tfr+ 


l/t-l 

+ ln(l — t) ln(^) 

V ^rnin / 

HO ~ T +)0 - r -)] 


V ~t 
1 - t 


In[(! -t/r+)(l -t/r_)], 


1 — t/r_ 
d'n = v Ln n=1...5. 

The formulae for C n and X/f are more compact than ([37]) and 


(45) 


C, = 

c 2 = 

C 3 = 

c 4 = 

C 5 = 
vf = 


—rjx + rjt + 5?/ — 2 xt + x + t + t 2 , 

\jx 2 — Ar] (6 — t 2 /r] + r / — 2x + 2 t), 

—5 + xt 2 /r] — 2rjx + ?/t + II77 — rf — 6 xt — x + 2a ; 2 — 3t + 5t 2 , 
—a:t 2 /?/ + 2r]x/t — 2 r]x — 10r]/t + 3r/t + 3?/ + 5rf /t 2 — 6 ?/ 2 /t 
—4a;t — 3x + 2a ; 2 + lit + 2t 2 , 

- \Jx 2 — Arj ^5 — 2t 2 jr\ — 5 r\/t — r] + 3tj , (46) 

Ci, 

5 


, U 1 ,2 ,1 1,0 3 . 5 9 

V 2o3 = ~-^+2 r lx/ t -'nx/t + -r]x--r]/t^-r]/t--7]t+-7i 

5 o lo 1 9 1 o 1 7 

^ A + 4*7 t + 1 ~xt + -xt + -x + -t 


}v 2 /t 2 


Vf = 


Vi 


V A = 


1 2 3 O 

+-t 2 - -t 3 , 

4 4 

\Jx 2 —Ar) (—6 + t 2 /i] — 7 ] + 2x — 2t) + 5 — xt 2 jr\ + 2r\x 
—7]t — 11?/ + 7] 2 + Qxt + x — 2x 2 + 3t — 5t 2 , 

\jx 2 — Ar\ (6 — t 2 /r] + 7 7 — 2x + 2t) + a;t 2 /?/ — 2r]x/t + 2?/a; 
-I - 1 077 /1 — 3 7]t — 3 T) — 5 ?/ 2 /t 2 + 6 ?/ 2 /t + 4a;t + 3a; — 2a ; 2 
— lit — 2 t 2 , 

\/x 2 -A7] + it 2 /?/ + i?//t + i?z - it) - ^oat 2 /?/ 
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1,1 ,1 11 3 2/ 9 5 

-~r]x/t + -rjx - rj/t - -rjt - —r) + -r/ /t + -xt + -x 

3 2 

+t - -t 2 . 


( 47 ) 


The coefficients T>^ are expressed by C n in the same way as D are expressed 
by C n in the set of equations i|iO|). 

The listed above results have been tested by comparison with the calcula¬ 
tion made earlier in simpler cases. One of the cross checks was arranged by 
fixing the mass of the produced lepton to zero. Our results are in this limit 
algebraically identical with those for the massless charged lepton|28|. |2iJ. On 
the other hand one can numerically integrate the calculated double differential 
distribution over x, with the limits given by the kinematical boundaries: 


fw++v/w+ dY dr. . 

/ T -77 ( x ,t;p,v) = -77 ■ 

> 2 y/fj ax at at 


(48) 


Obtained in such a way differential distribution of t agrees with recently 
published ||17| | analytical formula describing this distribution. This test is par¬ 
ticularly stringent because one requires two functions of three variables (t, p 
and rj) to be numerically equal for any values of the arguments. We remark, 
that for higher values of t only the region A contributes to the integral fl48[ ) 
and for lower values of t both the regions A and B contribute. This feature of 
the test is very helpful — the formulae for Fi(x,t ), which are different for the 
regions A and B can be checked separately. For completeness we quote the 
formulae derived in 0: 


dr 

d£ 


where 



-?)’ 
?o(t) - 


1 + 


2 a s 
3n 


m) 



2 a s 
3n 


•Pi (<) 


Mt) 

M*) 

m) 


4p 3 

(1 — p) 2 + 1( 1 + p) 

- 2t 2 

4p 3 

;(i-p) 2 -t(i+p); 

1 


+ A'lXw + 111 p + 

B, ^ + B 2 Y w + B-iYp + B4P3 111 p + B5P3 


(49) 


(50) 

(51) 

(52) 

(53) 
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and 


'b = 81n(2p 3 ) - 2Ini + [2Li 2 (w_) - 2Li 2 (tu+) + 4Li 2 (2p 3 /p + ) 

—4Yp ln(2p 3 /p + ) — Inp_ In w + + Inp + In w_] 2p 0 /p 3 , (54) 


A\ 

A. 2 

^3 

A4 

A 5 

Bi 

B 2 

B. i 

B 4 

b 5 


Fo(t), 

—8(1 — p) 1 + t — At 2 — p( 2 — t) + p 2 


-2 


3 + 6t — 21t 2 + 12t 3 - p( 1 + 12t + 5 1 2 ) 


+p 2 (ll + 2t) — p 3 


-6 


1 + 3t — At 2 — p{ 4 — t) + 3 p 2 


-2 [5 + 9t - 6t 2 - p(22 - 9t) + 5 p 2 

n it), 

-8(1 ~ p) [(1 ~ pf -t(l + p) 


-4(1 - p) A /t - 2(—1 + 3 p+ 15 p 2 - 5 p 3 ) 

+8(1 + p)t — 6(1 + p)t 2 , 

-4(1 - pf/t - 2(1 - p)(l - lip) + 6(1 + 3p)t, 
—6(1 — 3p)(3 — p) + 18t(l + p). 


(55) 


(56) 


The formula (|49| ) consists of two parts with transparent physical meaning 0 
The first term in the curly bracket containing the function J+ 0,1 (^) corresponds 
to the real t —> bW decay in the limit of vanishing b quark mass. The other 
one is related to a fictitious decay t —> bH + of a heavy quark t to the quark b 
and a charged Higgs particle H + . The formula j |32| | describing the first order 
QCD correction T\ (t) has been confirmed by many groups, cf. note added in 
[2D|| . On the other hand the formula for has been obtained only in 

and it still requires an independent check. Our calculation confirms the 


ref, 


33 


result (51) of ref. |33 


5 Moments of r energy distribution 

It has been argued in m 00 that moments of charged lepton energy 
distribution from 6-quark decay are valuable source of knowledge about the 
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physical parameters involved in the process and are only weakly affected by 
soft processes. Similar to the moments are the observables which are used in 
Id] to find the values of some nonperturbative parameters of HQET. 

Thus, in the notation of |T8| 


M n = 


El 


, dr 

lE r dE - 


r n = 


Mn 

M n 


(57) 


(58) 


where E min and E max are the lower and upper limits for r energy and M n 
involve both perturbative and nonperturbative QCD corrections to r energy 
spectrum. The nonperturbative 1/ml corrections to the charged lepton spec¬ 
trum from semileptonic B decays have been derived in the framework of 
IK)ET |[22|. [23] , 24] and even 1 /m\ corrections are known[35]. Up to order 


of 1 /ml the corrected heavy lepton energy spectrum can be written in the 
following way: 


1 dT 

Tn dx 




(59) 


where Ai and A 2 are the HQET parameters corresponding to the b quark 
kinetic energy and the energy of interaction of the b quark magnetic moment 
with chromomagnetic field produced by a light quark in the meson B. The 
functions can be easily extracted from the formula (2.11) in (2^. 


Following ref. |U| we expand the ratio r r 
2 a, 


= r,! 0 ) 


1 - 


37T 


■#>+4^’+4u 2> ). 

ml ml J 


where r ® is the lowest approximation of r r 


r (o) = 


m b 

2 




l+ri-p 


f 0 (x)x n dx 


Jv? p M x ) dx 


(60) 


(61) 


Each of the 6$ is expressed by integrals of the corresponding correction func¬ 
tion /W (x) and the tree level term fo(x) 


= 


/) 


1 +r)~P 

2v^ 


fd\x)x n dx 




x) dx 


Jity? P fo(x)x n dx 


12 $ P fo(x)dx ’ 


(62) 
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where the index i may denote any type of the described above corrections. 
The coefficients 8$ depend only on two ratios of masses of the charged lepton 
and the c quark to the mass m b . It is convenient to introduce the functional 
dependence of the 8$ on the masses in the following way: 


8^(m b ,m c ,m r ) = 8® 


f m b m c 
\m T ’ m b 


(63) 


Since the values of b and c quark masses are not known precisely, we calculated 
the 8’s in a reasonable range of them. We have chosen 4.4 GeV < m b < 
5.2 GeV and 0.25 < rn c /rn b < 0.35. In this range the dependence of the 
function 8$ r ^j on the variables can be approximated by a second order 
polynomial with the maximal relative error smaller then 1% . We propose: 


8(p,q) = a + b(p - po) + c(q - q 0 ) + d(p - p 0 ) 2 
+e(p - p 0 )(q - g 0 ) + f(q ~ go) 2 , 


(64) 


where p = m b /m T: p 0 = 4.75GeV/1.777GeV = 2.6730, g = m c /m b , g 0 = 0.28 
and the polynomial coefficients are fitted for each of the 8$ separately. Such 
a parameterization marks out the realistic masses of quarks: m b = 4.75 GeV 
and m c = 1.35 GeV, for which the coefficient 8$ = a$. 

The results of the fits of formula ( |64| ) to the calculated numerically co¬ 
efficients 8$(m b /m T ,m c /m b ) are listed in Table 1 for all the three types of 
8$ and n = 1... 5. With use of these results it is possible to estimate the 
order of relative correction to the ratio r n ([k]), including reasonable values 
of a s , Ai and A 2 . The value of A 2 may be easily determined from a mea¬ 
sured B - B* splitting, and a s should be of the order of a s (m b ). We as¬ 
sume a s ~ 0.3, A 2 = 0.12 GeV 2 , keep the mass of the b quark fixed to 
4.75 GeV and the mass of the c quark equal to 3.35 GeV. It was claimed 
in |Q that Ai may be constrained Ai = —0.35 ± 0.05 GeV 2 but we choose 
more conservative estimate —0.60 GeV 2 < Ai < —0.15 GeV 2 . Thus the per¬ 
turbative correction to r n /r ^ is about —0.0015 (—0.0075) for n — 1 (n — 5) 
whereas the nonperturbative corrections are larger: for the ’’kinetic energy” 
part 0.008 ± 0.004 (0.06 ± 0.03), n — 1 (n — 5) and for the ’’chromomag- 
netic” one —0.01 (—0.05), n = 1 (n — 5). Although the two nonperturbative 
terms partly cancel each other, the perturbative correction to the moments are 
smaller than the nonperturbative contributions. The ratios r n may be used to 
fix the value of Ai if suitable precise measurements are performed. 
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Table 1. Polynomial coefficients defined by the formula (^) characterizing 
how the three different types of corrections 5$ depend on quark masses. 


A. The perturbative correction. 



a 

b 

c 

d 

e 

/ 

IF 

0.0213 

0.0097 

-0.1059 

-0.0061 

-0.0202 

0.178 

~w 

0.0443 

0.0216 

-0.2213 

-0.0108 

-0.0454 

0.363 

~w 

0.0687 

0.0355 

-0.3457 

-0.0157 

-0.0808 

0.554 

~x(pT 

°4 

0.0946 

0.0512 

-0.4781 

-0.0193 

-0.1163 

0.746 

s(p) 

°5 

0.1219 

0.0684 

-0.6168 

-0.0259 

-0.1538 

0.934 


B. The nonperturbative correction corresponding to kinetic energy of the b 
quark. 



a 

b 

c 

d 

e 

/ 

F 

-0.500 

0.0 

0.0 

0.0 

0.0 

0.0 

IF 

-1.105 

-0.081 

0.220 

0.011 

0.016 

0.26 

T”" 

-1.823 

-0.254 

0.680 

0.035 

0.039 

0.84 

°4 

-2.664 

-0.526 

1.396 

0.074 

0.064 

1.77 

r(4) 

°5 

-3.633 

-0.904 

2.377 

0.131 

0.068 

3.11 


C. The nonperturbative correction corresponding to the interaction of the 
b quark spin with chromomagnetic field. 



a 

b 

c 

d 

e 

/ 

if" 

-1.731 

-0.732 

2.784 

0.212 

-0.498 

-1.71 

IF 

-3.565 

-1.561 

5.824 

0.437 

-0.975 

-3.45 

°3 

-5.496 

-2.473 

9.069 

0.681 

-1.506 

-5.13 

°4 

-7.517 

-3.453 

12.472 

0.961 

-2.141 

-6.72 

°5 

-9.621 

-4.485 

15.987 

1.279 

-2.955 

-8.21 
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6 Summary 


The results of this article can be summarized in the following way: 

(i) an analytical expression has been presented for the double differential dis¬ 
tribution dT /dx dt of the r lepton energy and the invariant mass of the leptons 
in b —> qrv T inclusive weak decay including first order perturbative QCD cor¬ 
rections. For m T = 0 the results of |28|, [29|] are rederived 

(ii) after numerical integration the distribution dY/dt is obtained in perfect 
agreement with the analytic result of [17] which has been derived in a com¬ 
pletely different way 

(iii) when combined with the non-perturbative QCD corrections calculated in 
the results of this article describe inclusive decays of beautiful 
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hadrons into r lepton and anything 

(iv) the energy spectrum and a few its moments are calculated for B into r 
transitions. 
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